Global and Local Information in Traffic Congestion 
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A generic network flow model of transport (of relevance to information transport as well as physical 
transport) is studied under two different control protocols. The first involves information concerning 
the global state of the network, the second only information about nodes' nearest neighbors. The 
global protocol allows for a larger external drive before jamming sets in, at the price of significant 
larger flow fluctuations. By triggering jams in neighboring nodes, the jamming perturbation grows 
as a pulsating core. This feature explains the different results for the two information protocols. 



The interplay between information networks and trans- 
port networks determine the flow on the latter. In par- 
ticular congestion formation, persistence and elimination 
depend on how information is spread across the trans- 
port system. In line with a number of recent studies 
III, 0, 0, m we model here these issues from a network 
theoretical perspective fd,^. Our aim is to study simple 
models to obtain qualitative as well as semi-quantitative 
insight into fundamental aspects of the dynamics of net- 
work congestion. Transport networks have indeed been 
object of many studies in recent years. In analogy to equi- 
librium fluid dynamics, steady-state descriptions [8] were 
proposed and studied in depth. Although they produce 
interesting results about local phenomena [1, [13, [HI, [lH 
they are not able to capture the inherent "intelligent" 
behaviour involved in route choice. Micro-level agent- 
based models [H, |l3, [HI solve the issue of individual 
choice, but require an enormous amount of detailed in- 
formations about the network (e.g. position and signaling 
of junctions, traffic lights synchronizing etc.). 
We propose here a minimal network flow model with 
unjammed/jammed nodes evolving under two different 
types of information distribution. In the first case we 
let each node have complete information {global model). 
In the second case each node only receive information 
concerning its neighborhood {local model). Nodes are 
thought of as stations or junctions in a traffic network, 
while links as roads or rail-tracks between two nodes. A 
node i is characterised by its threshold Ti = ftki, where 
ki is the node's degree and Q a positive constant, rep- 
resenting the maximum load it can support before jam- 
ming, the load Li representing the load (people, trains, 
cars..) present on the node, and the state Si, which is 
set to 1 or 0, respectively when the station is unjammed 
or jammed. 

When a node i jams, no incoming flow from neighbor- 
ing nodes is permitted until Li becomes smaller then Ti 
again. This mimicks the behavior of real systems, where 
a station can block due to too much incoming flow, ef- 
fectively cutting off its neighbors. Normal operation is 
restablished after sufficient off-loading. The equation for 
the state of node i is Si{t -\- 1) = & {Ti — Li{t)), where 
& is the Heaviside step function, and the outgoing flow 



Ff^^^ from node i to node j is 



min{Li{t),Ti) 
ki 



Sj{t){l-J{t)) (1) 



where J G [0, 1] represents the fraction of jammed 
stations that a station 'sees' or is informed of and 
the min function is used to put a superior limit on 
the outgoing flow from a station. Note that, for 
growing J, the outgoing flow from stations decreases 
in general. Also, the term Sj on the r.h.s. of ([T]) 
accounts for the impossibility of sending flow to a 
jammed station {Fij = if Sj = 0). The idea behind 
this choice is that a station, if informed that many 
stations are jammed, will try to gradually reduce its 
outgoing flux to avoid congesting the system any further. 

a. Initial Conditions We drive the system with two 
different mechanisms. In the first case, we introduce in 



tot 



toil 



the network a total load, L 
controlled by the parameter /3, that represents the filled 
fraction of the total network capacity. The total load is 
distributed in the network randomly, the only condition 
being avoiding to jam nodes from the beginning. Oper- 
atively, one can accomplish this by assigning to node i 
a load equal to a random fraction of the node's thresh- 
old and then iterating over the nodes until the load in- 
troduced in the network reaches the desired value. The 
initial conditions are thus: 



U{{)) = Xi 



S^{0) = 1 

0<Xi<Ti 



(2) 



where Xi a random variable G [0,Ti], so that "^Zi^i — 
PTtot- The system is then driven through random redis- 
tribution of load among nodes. 

In the second case, we place the entire load on one single 
node, denoted as the seed , and study how the system 
relax without any further drive. So the initial conditions 
are: 



Si{0) = l, Li{0)=0 yiy^seed 

Sseed — 0, Lseed — Ltot 



(3) 
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b. Dynamics The parameter J in ([T]) represents the 
fraction of jammed nodes that a station sees in the sys- 
tem. In the local dynamics, each station has informa- 
tion only about its nearest neighbors, thus J is node- 
dependent, J = Jl^'^t): 

^^W = r E (i-^.W) (4) 

where the sum is restricted to the neighborhood of node z, 
while in the global dynamics each station has information 
about the whole network, therefore producing an unique 
value for J for all the nodes at each time: 

J'^^it) = ^Y^il-Smit)) , (5) 

m 

We note that the completely jammed and unjammed 
states, Nj = N and Nj = 0, are absorbing states. In- 
deed, if the system reaches the state Nj = N we have 
jioc ^ J glob = 1 vi, implying Fij = Vz, j. Similarly, 
if Nj = 0, consider the load Li{to) on node i at time 
to. Since all the neighbors of i are unjammed, the to- 
tal outgoing flow from i is ^j^j. Fij = Li{to) while the 
maximum incoming flow is ^j^j. Fji = kia = Ti. Thus, 

Li{to + 1) = L{to) - E ^iJ + E ^ = 

jeh jeh 

and node i does not jam, since the condition for jamming 
is Li > Ti. 

We have seen that starting from Nj = the system can- 
not in itself trigger jamming perturbations. However, it 
is a very unstable state because, in presence of a jammed 
neighbor, a node keeps some fraction of its previous load 
on itself. This will move the node closer to the jam- 
ming threshold and thus making itself more susceptible 
to jamming. So if we produce in some way even a single 
jammed node, we expect the jamming perturbation to 
expand through a sort of chain reaction up to a station- 
ary state that depends on the amount of load introduced. 
This is a flrst indication that inhomogeneity in the sys- 
tem is the driving force behind jam propagation. 

c. Simulations We performed extensive simulations 
on random graph (RG) and scalefree networks (SF) with 
the number of nodes TV varying between 10^ and 5 x 10^ 
under both driving mechanisms. Realisations initialised 
as in ([2]) were driven by redistributing a fraction of a 
randomly chosen node's load to another randomly cho- 
sen node {Li ^ Li -\- cLj^Lj Lj{l — c)^c G (0,1)). 
When a jam appeared, the driving was suspended for 
the duration of the active phase {Nj ^ 0). Realizations 
initialized as in Eq. (|3]) relaxed to their stationary state 
without active driving. + Figure [T] shows the characteris- 
tic behaviors observed for the asymptotic jammed pop- 
ulation N j(t ^ oo). For < 0.75 the system does not 
show a stationary population of jammed nodes For 
higher values, we observe N j ^ for different values 
of depending on the topology and dynamics, but, in 




FIG. 1: NjYS (3 plot for global and local models on SF (full 
line) and RG (dashed line) networks {N = 1000) with seeded 
and uniform (inset) initial conditions. Error bars are stan- 
dard deviations of fluctuations of Nj around its mean value, 
averaged over the realizations at a fixed f3. Data shown are 
obtained over 100 realizations of the RG and SF networks. 



striking contrast to the naive expectation Nj ^ A/", as 
P ^ 1 Nj remains in he vicinity of N/2 for a broad 
range of (3 values. Indeed the applied load is suflftcient 
to jam the entire network. However the load is trapped 
in the jammed nodes and being redirected from there. 
This prevents a uniform redistribution of the load across 
the network keeping Nj far below its maximal possible 
value. We notice that in the local model the increase in 
Nj sets in for lower values of P than in the global model. 
Also, the SF networks appear to be more susceptible of 
developing a stationary jammed population than the RG 
network. There reason for this is the higher abundance 
of low degree nodes in SF networks than in RG networks. 
The low degree nodes are more easily jammed and pro- 
vide the core needed for seeding a larger jamming pertur- 
bation. Indeed, Nj{t) under the seeded driving (flgE]) in- 
creases through a sequence of jumps and quasi-plateaus, 
these are the signature of an oscillating core mechanism 
(flgE]): the initial jammed seed makes the neighboring 
nodes more vulnerable to jamming by sending flow to- 
ward them, while not accepting flow from them. This 
brings the seed's neighbors closer to their threshold. If 
there is enough load on the seed node, it will eventu- 
ally jam its own neighborhood (flg. [3ti); the neighbors 
will relax onto their unjammed neighbors, but this will 
make the second neighbors vulnerable or even jam (flg. 
[3t)) ; this mechanism can continue as long as the core has 
flow to send. The growth of the core can end in two 
scenarios, a) the seed node's neighbors manage to dis- 
tribute the load without jamming their own unjammed 
neighbors and so the load is distributed freely outward 
(flg. [3ll), or b) the second neighbors jam and while they 
relax, they jam the seed's flrst neighbors again thus ere- 
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ating a bigger, more stable jammed core (figl3l:). The 
quasi-plateaus in Nj{t) correspond to periods of growth 
of the vulnerable population, and are followed by sharp 
jumps, where the jammed core grows by rapidly invad- 
ing the vulnerable nodes. Another interesting feature 
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FIG. 2: Time series plot for global (green) and local (black) 
models in a RG network with /3 = 0.9, averaged over realiza- 
tions. 

is that the local model exhibits much narrower fluctu- 
ations {ai ^ 30 for N = 10^) than the global model 
{ag 150). These fluctuations are another signature of 
the core expansion mechanism: when the systems reaches 
its stationary Nj value, it keeps fluctuating around the 
value due to the constant overshooting between jammed 
and unjammed nodes. The more violent fluctuations ob- 
served under the global dynamics are indeed consistent 
with the core mechanism: the global dynamics, through 
its global flow suppression, produces smaller "packets" 
of flow, that allow nodes to stay nearer to their thresh- 
olds than under the local dynamics, and therefore pro- 
duces a larger population of nodes susceptible of jam- 
ming, thereby generating bigger fluctuations. For this 
same reason, the global dynamics develops a stationary 
jammed population for /3 values larger than for the local 
dynamics but exhibits also a steeper increase with (3. 

d. Fokker- Planck Rate Model Our deterministic 
flow model is computationally heavy to simulate and an- 
alytically hard to describe. Although the model is deter- 
ministic, the fluctuations due to spatial inhomogeneity 
suggest it may be possible to understand observed be- 
havior by use of a set of simple transition processes. We 
separate the population of nodes in to three types: un- 
jammed (U), vulnerable (V) and jammed (J), with rate 
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FIG. 3: Jammed component expansion mechanism. 




FIG. 4: An example of fluctuations PDF for RG for [3 
0.91{N = 1000). 



equations, 

/7(n,t) + J(m,t) A F(n,t) + J(m,t) (7) 

V{n,t) ^ J{m,t) A J{n,t)^J{m,t) (8) 

J{n,t)^U{m,t) A U{n,t)^U{m,t) (9) 

Instead of treating individual nodes, we consider degree 
classes and introduce Uk,Vk,Jk as the probabilities of 
choosing randomly a member of a /c-degree node belong- 
ing to one of the three types. The rate terms are very 
intuitive and can be written down directly from the rate 
equations, considering the contributions from the differ- 
ent degree classes. For example, the contribution of (|9]) 
to Jk is of the form Ti{k,t) = ^Efc^=i P{Vk{t), Jk'{t)). 
Later on, we will set v = 1 as deflnition of vulnerable 
population. For the moment we neglect degree-degree 

1000^ ^ ^ ^ ^ ^ ^ 1 

□ 
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FIG. 5: SF simulated data {dots) and fit {line) obtained using 
eq. [14] for local (green and blue) and global (violet and red) 
models. 

correlations and write the probabilities P(a, h) as the 
probability of randomly choosing a node belonging to 
a and a second node belonging to h times the proba- 
bility of having a link between the two. For example, 
PiVk^Jk') = Hj (AT-i) 7 other words the process 
rate times the number of relevant couples VJ. If we 
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now substitute the F terms inside the FP equations for 
the populations we obtain, for each degree class /c, the 
set of equations: 



N{N -1)2M 



(10) 



k' = l 

Uk{t + 1) = Uk{t) + 

kmax 



k' = l 

Vk{t + 1) = Vk{t) + 

kmax 



k' [vVk{t)Jk'{t) - lJk{t)Uu'{t)] 
k 

N{N -l)2M ^ 
k' hJk{t)Uk'(t) - aUk{t)Jk'{t)] 
k 



(11) 



N{N-1)2M 
J2 k' [aUk{t)Jk'{t) - uVk{t)Jk'{t)] 



(12) 



k' = l 



with Jj, ^Uk ^Vk = Pk and ^j.Pk = 1- Performing the 
sum over k'we obtain the average degree for the respec- 
tive populations. The (|lQtilip cannot be solved easily. 
We restrict to the case (kj) = (ku) = {ky), a strong as- 
sumption that is however confirmed by the simulations. 
So for the stationary state we find: 



Jk = -Vk Uk = -Vk yfe = — ^ ^ 
7 a l+a+^ 



(13) 



Assuming a simple proportionality between /3 and the 
"jamming rate" a ex {P — Po) , we can compare the 
Fokker-Planck prediction 



iP - /3o) 



7+(7+l)(/9-/9o) 



(14) 



with the simulated data. Figure [5] shows this compari- 
son and we can easily see that the Fokker-Planck solu- 
tion reproduces well the behavior of the simulated data 



up to :^ 1. In particular, the fit correctly identifies 
the value of Po and reproduces the steep transition from 
no stationary non-zero Nj to a significant non-zero Nj. 
As expected, eq. [TH does not predict the second jump 
in Nj N: the approximation of uncorrelated degrees 
and constant rates break down. Indeed, eqs. ([HE]) are 
not valid anymore as the system is better described by 
only two competing population V and J. Moreover, our 
mean-field-like rate model does not account for random 
fiuctuations that could bring the system in the absorbing 
Nj = N state. 

e. In conclusion we have identified how a jammed 
core is able to push neighboring nodes to their thresh- 
olds and eventually invade them. We also demonstrated 
that the network topology and the different information 
regimes (local and global) play a significant role only at 
the onset of jams (by shifting threshold load level /3), 
but do not significantly infiuence the rest of the evolu- 
tion, which is qualitatively well described by a model of 
3 competing populations, mimicking the unjammed (U), 
vulnerable (V) and jammed (J) populations of the fiow 
model. 

Of relevance to traffic control we have shown that un- 
der global dynamics the jam picks up at a larger driving 
load than when only local information is used. However, 
in our model global traffic management exhibits much 
larger ffuctuations and a steeper increase in number of 
jammed node. Indeed, for real systems, our model sug- 
gests that under heavy loads it is better to let the sys- 
tem evolve under local information rather than attempt- 
ing global management, since jams sweep a much larger 
part of the system in the latter case. Moreover, since 
congested regions enhance the jamming perturbation by 
making their neighbors more susceptible to jamming, it 
is better to reduce the number of links connecting the 
congested regions to non-congested ones (in some sense, 
quarantine them) rather than letting them communicate 
with the neighboring regions and thereby risk the trig- 
gering of new jams. 
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